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We propose an explanation of the mid-infrared peak observed in the optical conductivity of the 
Bechgaard salt (TMTSF^PFo in terms of electronic excitations. It is based on a numerical calcu- 
lation of the conductivity of the quarter-filled, dimerized Hubbard model. The main result is that, 
even for intermediate values of U/t for which the charge gap is known to be very small, the first 
peak, and at the same time the main structure, of the optical conductivity is at an energy of the 
order of the dimerization gap, like in the infinite U case. This surprising effect is a consequence of 
the optical selection rules. 

PACS numbers: 71.10.Pm, 7I.27.+a, 78.20.Bh, 78.30.Jw 
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I. INTRODUCTION 

It is by now well established experimentally that there 
is a well defined structure in the mid-infrared optical con- 
ductivity of the Bechgaard salt (TMTSF) 2 PF 6 §-§. Ac- 
cording to the most recent data of Dressel et al [3|, this 
structure moves from about 1000 cm -1 at room temper- 
ature to 200 cm -1 at 20 K, while its intensity increases 
upon lowering the temperature. Given the width and the 
intensity of this structure, it seems difficult to explain it 
in terms of phonons, and the most natural thing to do it 
to look for an explanation in terms of electronic excita- 
tions. 

It is not so easy however to understand why there 
should be an absorption due to electronic transitions 
in that energy range. To see that, let us consider the 
simplest description of the Bechgaard salt that contains 
the essential physics, namely the quarter-filled, dimerized 
Hubbard model described by the following Hamiltonian: 

H=-t x + 

i even,(T 

-*2 ( C ]+la C «r + h - c -) + U Y1 n ^ n H (!) 

i odd. a - i 

where ct , c ia create and annihilate holes in the HOMO 
of the TMTSF molecules. For (TMTSF) 2 PF 6 , reason- 
able parameters are ii=250 meV, tz/ti = -9 and U/t\ = 
5. Note that a Luttinger liquid description would not be 
of much use here given the energy range we are interested 
in. 

Now, let us see which electronic transitions could be 
present in such a model. If the interactions are ignored, 
the only allowed transitions are the vertical transitions 
from the lower band to the upper b and. T hey give rise 
to a continuum starting at Ao = %\Jt\ +t\ ~ 670 meV, 
which is roughly one order of magnitude too large. In- 
cluding U and assuming that U jt\ = 5 will induce tran- 
sitions at energy U ~ 1.25 eV, which is even larger. In 



fact there is only one characteristic energy scale that has 
the right order of magnitude, namely the dimerization 
gap Ad = 2(ti - t 2 ) = 50 meV. This led Pedron et al 
Pj to assume that the double occupancy of a site must 
be excluded, or equivalently that U is infinite. Then the 
charge carriers can be described as a half-filled system 
of spinless fermions, in which case the vertical transition 
has an energy of Ad- There is a problem however with 
that explanation. If U is infinite, then we know that 
the charge gap A p is equal to the dimerization gap Ad 
But such a large charge gap is inconsistent with the 
experimental fact that the conductivity remains metal- 
lic down to very low temperatures. In fact, the metallic 
character of the conductivity puts an upper bound on the 
charge gap which is consistent with the above-mentioned 
parameters j^]. 

In this paper, we show that the structure observed in 
the optical conductivity is actually consistent with the 
model of Eq. (1) with <i=250 meV, t 2 /h = .9 and 
U/t-y = 5. The main point is that the charge excita- 
tions with energy around A p responsible for the metallic 
conductivity cannot be seen in the optical conductivity 
because of the optical selection rules. The first allowed 
transition occurs at an energy of the order of the dimer- 
ization gap Ad, which is much larger than A p for these 
values of the parameters. 

The paper is organized as follows. In section II, we first 
discuss the general properties of the optical conductivity 
in the model of Eq. (1). In section III, we describe the 
numerical method and we derive the form of the finite- 
size corrections we can expect for the position of the first 
peak of the optical conductivity. In sections IV and V, 
we present the results for large and small dimerization 
respectively. Finally, a discussion of the results in con- 
nection to (TMTSF) 2 PF 6 is given in section VI. 
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II. GENERAL PROPERTIES OF THE 
CONDUCTIVITY 

Let us start with a qualitative description of the optical 
conductivity in various cases of the Hamiltonian of Eq. 
(1) at quarter-filling. They are illustrated in Fig. ||. 

a) Non interacting electrons on a non-dimerized lattice 
(U = 0, t<x = t\): In that trivial limit, the only contribu- 
tion is of course the zero frequency Drude peak. 

b) Interacting electrons on a non-dimerized lattice 
(U 5^ 0. t2 = t\): There is still a Drude peak because the 
system remains metallic, but there is also some incoher- 
ent spectral weight around u) = U. 

c) Non interacting electrons on a dimerized lattice 
(U = 0. t% ^ t\): The dimerization splits the disper- 
sion into 2 bands with the equation: 



E{k) = ± 



f? + *| 



2t\ti cos(2fcet) 



(2) 



where 2a is the lattice parameter. The conductivity ex- 
hibits a Drude peak since the system is metallic, and 
a continuum between Ao and 2(*i + t%) due to vertical 
interband transitions. See Fig. |[ 

d) Interacting electrons on a dimerized lattice (U =/= Q, 
t-2 i= ti)- In that case the repulsion opens a gap because 
the lower band is effectively half-filled. So the system 
becomes insulating and the Drude peak disappears. Now 
the incoherent and Drude weight are related by the sum- 
rule 



a(uj)duj = < — T > 

n 2L 



(3) 



where < T > /L is the expectation value of the kinetic 
energy per-site in the ground state. This quantity is not 
dramatically reduced by the interaction j7j, and the lost 
Drude weight must be redistributed as an incoherent, low 
energy background. Part of the weight can also migrate 
to the upper Hubbard band at an energy of order U. 

If U is very large, the low energy, incoherent weight 
must be located around Ajj to become the zone- 
boundary interband transition of spinless fermions in the 
U = +oo limit described by Pedron et al [||. In the fol- 
lowing our main goal is to determine whether for realistic 
parameters, i.e. for intermediate values of U/t\, there is 
still a dominant low energy structure in the conductivity, 
and where it is located. 



In this formula, Eq is the ground state energy, n labels the 
excited states and j is the paramagnetic current operator. 
The weight of the Drude peak D is calculated from the 
Kohn relation §-|ll): 



(4) 



where Eo(<j)) is the ground state energy as a function 
of the twist in the boundary conditions, while the inco- 
herent part can be obtained as a continued fraction using 
Lanczos algorithm || . Extensive calculations of that sort 
have been performed for 2D models in the context of high 
T c superconductors f|. 

We have performed exact diagonalizations on finite size 
clusters with L = 8, 12, 16, 20 sites. To reach the ther- 
modynamic limit (L = +oo) we have tried to perform a 
finite-size scaling of the results. The finite size correc- 
tions turn out to be very large concerning the location f2 
of the first peak in the conductivity for the following rea- 
sons: In the infinite U limit, the charge and spin variables 
are decoupled (see Ogata and Shiba Jl2|). The charge 
part is that of spinless fermions with twisted boundary 
conditions and this twist is given by the momentum of 
the spin part. Now, for U = +oo there is no energy 
associated to spin excitations, and the spinless fermions 
are free to choose the boundary conditions that minimize 
this energy, namely antiperiodic boundary conditions. In 
that case, the allowed values of the momentum are given 
by k — (2v + l)n/La (see Fig. ||). The zone boundary 
k = ti /2a does not belong to the allowed momenta and 
the larger occupied k is at it /2a — tt/ La. So the first peak 
in the conductivity is located at 



Q = A D (L) = ytf + f| + 2ii*2 cos(tt - 2n/L) (5) 

In other words, there is a finite-size correction to the 
known result O = Ao in the thermodynamic limit. In 
the case of large but finite U, the spin energy is still 
negligible compared to that of the charge degrees of free- 
dom, and the charge part keeps the same boundary con- 
ditions. This suggests that the finite-size corrections to 
Q? will scale according to cos(7r — 2tt/L). We shall see 
that this scaling form is actually remarkably accurate 
down to rather small values of U / 1\ . 



IV. LARGE DIMERIZATION 



III. THE NUMERICAL METHOD 

Our aim is to calculate the optical conductivity of the 
model at zero temperature. Starting from its definition 
as the current-current correlation function, the real part 
of the conductivity can also be written H : 



a(w) = DS(lj) 



E n — Eq 



-8{u-E n +E ) 



For clarity, let us start with a case where all the impor- 
tant features of the conductivity can be best seen in spite 
of the limitations due to the finite size of the clusters. In 
units of *i we choose t 2 = 0.5 and study various repul- 
sions. On Fig. | we show the influence of the cluster size 
on the incoherent spectum for a given repulsion. This 
incoherent part exhibits a sharp peak which is shifted 
toward low frequencies when the size increases, while its 
relative weight is more or less constant. Figure shows 
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the effect of the repulsion {U/ti change from to 10) for 
L=16. The main peak increases with size and approaches 
Ad(L) for large U as expected. 

Because of the dimerization, the lower band is half 
filled and the repulsion induces an insulating behavior in 
the thermodynamic limit. So the Drude weight given by 
Eq. (5), which does not vanish for finite systems, should 
scale to zero upon increasing the size of the cluster. Our 
numerical results are effectively consistent with a vanish- 
ing Drude weight in the thermodynamic limit, although 
the finite size scaling function is not very clear (see Fig 

I). 

We have also checked the sum-rule of Eq. (3). The fi- 
nite range of integration could bring problems in practical 
computations, but in the case of quarter-filled systems, 
there is essentially no weight at high frequencies [ p^|Jl4| 
and one can safely stop the integration at to = 2U. In all 
cases we found that the sum-rule was satisfied with an 
accuracy better than 1%. 

In order to show that the incoherent conductivity ex- 
hibits a well defined structure at low energy, we have 
calculated the relative weight of the first peak with re- 
spect to the total incoherent part. The results are given 
in Table [j]. This proportion is insensitive to the size and 
increases with the repulsion. This suggests that the first 
peak will dominate the incoherent part in the thermody- 
namic limit, even for intermediate repulsions. 

To find the location of this peak, we have tried sev- 
eral finite-size scaling. The only way to obtain a good 
scaling is to plot f2 2 versus cos(7r — 2ir/L). Some results 
are presented in Fig. [| The extrapolated values are 
given in Table II. In the present case, the location of the 
first peak roughly follows the charge gap A p , which is 
relatively large for this value of the dimerization. 

V. SMALL DIMERIZATION 

We now turn to a smaller dimerization, namely t% — 
0.9ti, having in mind intermediate repulsions U ~ 5t\. 
The incoherent conductivity <Ji nc calculated for L=20 is 
shown in figure [?]. As in the previous case Ui nc is still 
dominated by its first peak. The relative importance of 
this peak increases from 34.3% on 8 sites to 53.1% on 20 
sites. 

In the present case <Ji nc represents a rather small part 
of the oscillator strength (from 2 % for 8 sites to 4.5 % for 
20 sites). Because of the weak dimerization, the charge 
gap is small (A p = 0.02ti), and the corresponding length 
vf/ 'A p is much larger than the sizes we can reach. So the 
Drude peak takes almost all the oscillator strength in our 
simulation. But this peak will certainly disappear in the 
thermodynamic limit since the system is insulating. So 
what really matters is whether the relative weight of the 
first peak be sizable, which it is. 

Concerning the location of the first peak, we found 
that scaling fl 2 with cos(7r — 2ir/L) was still very accu- 



rate and for U — 5ti this scaling leads fl/ti = 0.17 in the 
thermodynamic limit. Now, for t% = -9ii and U = 5t±, 
the charge gap A p equals 0.02 and the dimerization gap 
Ad equals 0.2. So we found that the first peak appears 
at an energy which is of the order of the magnitude of the 
dimerization gap, and that there is no weight at energies 
corresponding to the charge gap. How can we understand 
this result? If the charge gap is identified to the lowest 
excitation of the system which leaves the total spin S 
unchanged, the difference between and A p should be 
directly observable on the excitation spectrum. Such a 
spectrum is plotted in Fig. |^. The lowest charge exci- 
tation occurs at the border of the Brillouin zone, while 
f2(L) corresponds to a vertical excitation keeping S=0. 
Besides, the first singlet state coupled to the groundstate 
by the current operator is not the lowest S=0 excited 
state but the third one jl(J. This is presumably due to 
the fact that the current operator being odd under the 
inversion, only odd states can be coupled to the ground- 
state, which is even. Q(L) and A p are thus clearly dif- 
ferent for L finite. Besides, they have very different scal- 
ings: On one hand A„ goes to a very small value like 
l/L as shown on Fig g. The exact value 0.02 provided 
by a perturbative calculation can actually not be identi- 
fied by such a scaling. On the other hand, f2(L) 2 scales 
as cos(7r — 2tt/L) toward 0.17. 

So our results indicate that the first peak of the opti- 
cal conductivity is located at an energy much larger than 
the charge gap. We cannot exclude that, upon increasing 
the size, other states will appear that have the right sym- 
metry to be coupled to the groundstate by the current 
operator. Such states could for instance involve several 
excitations of momentum n and of energy of order A p 
leading to a small threshold at 2A p in the optical con- 
ductivity. But according to the increase with L of the 
relative weight of the first peak (from 34.3% for 8 sites 
to 53.1% for 20 sites), the feature at fl ~ 0.17ii should in 
any case remain the prominent structure of the conduc- 
tivity in the thermodynamic limit. Finally, the difference 
with the case treated in the previous section presumably 
comes from the fact that the excitations have very lit- 
tle dispersion when the dimerization is large because the 
bands are flat. 

So it seems that the explanation proposed by Pedron et 
al |H is essentially valid: For intermediate values of U /t-y 
and small dimerization, the conductivity has a large peak 
around the dimerization gap Ad although the charge gap 
is very small. 

VI. DISCUSSION 

In conclusion, we have shown that the model of Eq. (1) 
with reasonable parameters leads to a peak in the opti- 
cal conductivity at an energy of about 0.17ii ~ 40 meV. 
This energy has the right order of magnitude. So this 
provides a good candidate to explain the mid-infrared 
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structure observed in the Bechgaard salt (TMTSF) 2 PF 6 . 
What about the fine details of this mid-infrared struc- 
ture? As stated in the Introduction, this peak moves to- 
ward lower energies when the temperature is decreased, 
and its weight increases. Both features are actually quite 
natural consequences of the present explanation. 

Let us first consider the fact that the peak moves to 
lower energies when the temperature is lowered. The 
Hamiltonian of Eq. (1) , which is purely ID, is not an ac- 
curate description of the electronic structure of the Bech- 
gaard salts at low temperatures where 2D effects start 
to dominate. Looking at the actual 2D band structure 
calculated on the basis of the structure deduced from 
X-Ray measurements performed at 300 K and 4 K, one 
can see that there is a clear evolution leading to a much 
smaller dimerization gap at low temperatures. Whether 
this gap still yields a structure in the optical conductiv- 
ity when the system must be considered as 2D cannot 
be concluded on the basis of the present calculation and 
should be checked independently, but it seems likely that 
the structure we have calculated will move smoothly to 
lower energies. 

Concerning the intensity, the data of Dressel et al ||] 
suggest that the increase of weight of that structure upon 
lowering the temperature is accompanied by a decrease of 
the weight of the background, which can be interpreted 
as a broadened Drude peak. Although the model we are 
looking at is insulating at zero temperature, we expect a 
Drude peak to be restored at temperatures comparable 
to the charge gap A p ~ .02^ ~ 40 K. It seems plausi- 
ble that the weight of this peak will thus increase with 
temperature in the range 50-100 K. Note that the fact 
that the system is metallic below 40K is not incompati- 
ble with A p ~ 40 K because the physics is certainly 2D 
below that temperature |l5|], and the ID description does 
not apply anymore. It would be interesting to check if 
this evolution can be confirmed by Quantum Monte Carlo 
simulations. 

We acknowledge very useful discussions with L. De- 
giorgi, M. Dressel and X. Zotos. We are especially grate- 
ful to D. Poilblanc for very useful comments on several 
aspects of the calculation of the conductivity by exact 
diagonalization of finite clusters. The numerical calcu- 
lations have been performed on the C-94 and C-98 of 
IDRIS (Orsay). 
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TABLE II. Comparison between the charge gap and the 
location of the first peak in the optical conductivity for large 
dimerization (£2 = 0.5ii) and various repulsions. The value 
of A p are taken from Ref. 8. 



FIG. 1. Schematic picture showing the influence of the re- 
pulsion and of the dimerization on the str ucture o f the con- 
ductivity. A D = 2(ti - ti) and A = 2yjt\ + t\. Fig. (d) 
is typical of large values of U. For clarity, the 5 peaks of the 
incoherent part of the conductivity have been broadened. 

FIG. 2. Dispersion in the non-interacting case with dimer- 
ization. (1) Vertical transition in the non interacting case 
at uj — Aq. (2) Finite size dimerization gap Ad(L — 16). 
(3) Dimerization gap. The circles (resp. triangles) are the al- 
lowed wave vectors for L=16 and periodic (resp. antiperiodic) 
boundary conditions. 

FIG. 3. Evolution of the optical conductivity with the size 
of the cluster, for £2 = 0.5£i, U = 5£i and e = 0.05£i (width 
of the S function). 
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FIG. 4. Evolution of the optical conductivity with the on 
site repulsion U for ti = 0.5£i, L = 16 and e = 0.05£i (width 
of the S function). 



FIG. 5. Finite size scaling of the relative weight of the 
Drude peak (D(L)/(Tve 2 < —T > /£)) for £2 = 0.5£i. 
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FIG. 6. Finite size scaling of the square of the frequency 
Q(L) (in units of *i) of the first peak of the incoherent con- 
ductivity. 

FIG. 7. The optical conductivity for * 2 = 0.9<i, U = 5*i, 
L = 20 and e = 0.05*i (width of the 8 function). 



FIG. 8. The first excitations versus the impulsion for 
*2 — 0.9, U = 5, L = 20. The impulsion k is related to v 
by k = —it/L + 2-kv/L because of antiperiodic boundary con- 
ditions. The circles (resp. triangles) represent states with a 
total spin S=0 (resp. S / 0). (1) First excitation coupled 
to the ground state by the current operator. (2) First singlet 
excitation. 



FIG. 9. Finite size scaling of the lowest charge excitation 
for ti = 0.9 and U = 5. 
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